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I. INTRODUCTION 

The search for CP violation (CPV) in charm decays 
provides a sensitive probe of physics beyond the Standard 
Model (SM). Owing to its suppression within the SM, a 
significant observation of direct CPV in charm decays 
would indicate the possible presence of new physics (NP) 
effects in the decay processes. In a previous article [T], we 
reported a precise measurement of the CP asymmetry in 
the D ± — > mode, where the measured asymmetry 

was found to be consistent with the value expected from 
indirect CPV in the K° system. 

The LHCb and CDF Collaborations have recently re- 
ported evidence for CPV in charm decays by measuring 
the difference of CP asymmetries in the D° — > K + K~ 
and D° — > 7r + 7r~ channels [21 [3], which is mainly sensitive 
to direct CPV. The size of the world average direct CP 
asymmetry difference, (—6.56 ±1.54) x 10 -3 [J], suggests 
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either a significant enhancement of SM penguin ampli- 
tudes or of NP amplitudes (or both) in charm decays [5] . 
Improved measurements of the CP asymmetries in the 
individual two-body modes, along with measurements in 
other channels, are needed to determine the nature of the 
contributing amplitudes. 

We present herein measurements of the decay rate CP 
asymmetry, Aqp, defined as 



A C p = 



T(D+ 



0) 



f)-r(D 



/) 



T(D + 



/) + r(i> 



7)' 



(i) 



in the decay modes D ± — > K®K ± , Df 
and Df -+ K^ir . Previous measurements of Aqp in 
these channels have been reported by the CLEO-c [S] 
and Belle Collaborations [7]. As for the Aqp measure- 
ment in D ± — > K^^, we expect an Acp asymmetry of 
(±0.332 ± 0.006)% [5] resulting from CPV in K° - K° 
mixing [9 . The sign of the -fT°-induced asymmetry is 
positive (negative) if a K° (K°) is present in the cor- 
responding tree-level Feynman diagram. Because it is 
identified by its 7r + 7r~ decay, the intermediate state is 
a coherent mix of and amplitudes. It has been 
shown in Ref. 10 j that the K® — K® interference term 
gives rise to a measured CP asymmetry that depends on 
the range in proper time over which the decay rates are 
integrated, and on the efficiency for the reconstruction of 
the intermediate state as a function of its proper flight 
time. 
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For this analysis we employ a technique similar to that 
used for our measurement of CPV in the D ± — »• K^tt^ 1 
mode pp. As a result, reference to our previous publica- 
tion is given for the description of some of the analysis 
details. 



II. THE BABL4R DETECTOR AND EVENT 
SELECTION 

The data used for these measurements were recorded 
at or near the T(4S) resonance by the BaBAR detector 
at the PEP-II storage rings, and correspond to an inte- 
grated luminosity of 469 fb _1 . Charged particles are de- 
tected, and their momenta measured, by a combination 
of a silicon vertex tracker, consisting of 5 layers of double- 
sided detectors, and a 40-layer central drift chamber, 
both operating in a 1.5 T axial magnetic field. Charged- 
particle identification is provided by specific ionization 
energy loss measurements in the tracking system, and 
by the measured Cherenkov angle from an internally 
reflecting ring-imaging Cherenkov detector covering the 
central region of the detector. Electrons are detected 
by a CsI(Tl) electromagnetic calorimeter. The BABAR 
detector, and the coordinate system used throughout, 
are described in detail in Refs. [TTJ [15]. We validate 
the analysis procedure using Monte Carlo (MC) simu- 
lation based on Geant4 [13]. The MC samples include 
e + e _ — > qq (q = u,d,s,c) events, simulated with JET- 
SET [13] and BB decays simulated with the EvtGen gen- 
erator 15j. To avoid potential bias in the measurements 
we finalize the event selection for each channel, as well as 
the procedures for efficiency correction, fitting, and the 
determination of the systematic uncertainties and possi- 
ble biases in the measurements, prior to extracting the 
value of Aqp from the data. 

Signal candidates are reconstructed by combining a 
candidate, reconstructed in the decay mode — > 7r + 7r~ , 
with a charged pion or kaon candidate. A K® candi- 
date is reconstructed from two oppositely charged tracks 
with an invariant mass within a ±10 MeV/c 2 interval cen- 
tered on the nominal K® mass [5] , which is approximately 
±2.5 a in the measured mass resolution. The x 2 - 
probability of the ir + ir~ vertex fit must be greater than 
0.1%. Motivated by MC studies, we require the mea- 
sured flight length of the candidate to be at least 
three times greater than its uncertainty, to reduce com- 
binatorial background. A reconstructed charged-particle 
track that has px > 400 MeV/c is selected as a pion or 
kaon candidate, where pt is the magnitude of the mo- 
mentum in the plane perpendicular to the z axis (trans- 
verse plane). In our measurement, we require that a pion 
candidate not be identified as a kaon, a proton, or an elec- 
tron, and that a kaon candidate be identified as a kaon, 
and not as a pion, a proton, or an electron. Identifica- 
tion efficiencies and misidentification rates for electron, 
pions, kaons, and protons with 2 GeV/c momentum in 
the laboratory frame are reported in Table [I] The crite- 



TABLE I: Identification efficiencies and misidentification 
rates for electron, pions, kaons, and protons with 2 GeV/c 
momentum in the laboratory frame. The values for kaons 
on the third row refers to the identification criterion used to 
reject kaons from the pion sample, while the values on the 
fourth row to the criterion used in the kaon selection. 



Particle 


Efficiency [%] 


Misid. rate [%] 




A± 


e± 


91 


0.04 


< 0.2 




88 




1 


A' (applied to tt ) 


91 


1 




A± (applied to AT*) 


99 


8 




P ± 


80 


0.2 


0.2 



ria used to select pion or kaon candidates are very ef- 
fective in reducing the charge asymmetry from track re- 
construction and identification, as inferred from studying 
the data control samples described below. A vertex fit to 
the whole decay chain, constraining the D,^ production 

vertex to be within the e + e~ interaction region, is then 
performed [IB]. We retain only candidates having 

a x 2 -probability for this fit greater than 0.1%, and an 
invariant mass m(K®h), h — ir,K, within a ±65MeV/c 2 
interval centered on the nominal Z)^\ mass [5], which is 

approximately equivalent to ±8 a in the measured 
mass resolution. 

We require further that the magnitude of the Df can- 
didate momentum in the e + e~ center-of-mass (CM) sys- 
tem, p* , be between 2.6 and 5.0 GeV/c, in order to sup- 
press combinatorial background from BB events. For the 
D ± — > K t ^K ± mode, the MC simulated sample shows 
that retaining candidates with p* between 2.0 and 5.0 
GeV/c allows signal candidates from B- meson decays, 
without introducing an excessive amount of combinato- 
rial background. Assuming that CPT is conserved, there 
is no contribution to Aqp from CP violation in B me- 
son decays from Standard Model processes. Additional 
background rejection is obtained by requiring that the 
impact parameter of the candidate with respect to 
the beam-spot |llj . projected onto the transverse plane, 
be less than 0.3 cm, and that the proper decay time, 
t xy , be between —15 and 35 ps. The decay time is mea- 
sured using L xy , defined as the distance of the DjK decay 
vertex from the beam-spot projected onto the transverse 
plane. 

In order to further optimize the sensitivity of the Acp 
measurements, we construct a multivariate algorithm, 
based on seven discriminating variables for each D, *, 
candidate: t xy , L xy , p* , the momentum magnitude and 
component in the transverse plane for the K® candi- 
date, and also for the pion or kaon candidate. For the 
D ± -> KgK^ 1 and Df -» K^K ± modes the multivariate 
algorithm with the best performance is a Boosted De- 
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cision Tree Q7], while for the Df -> K^ir* mode the 
best algorithm is a Projective Likelihood method [17]. 
The final selection criteria, based on the outputs of the 
multivariate selectors, are optimized using truth-matched 
signal and background candidates from the MC sample. 
For the optimization, we maximize the S/y/ S + B ratio, 
where S and B are the numbers of signal and background 
candidates with invariant mass within ±30 MeV/c 2 of the 
nominal Z)^\ mass, which is approximately ±3 a in the 
measured mass resolution. 



asymmetry (A FB ), and a detector- induced component. 
We measure A FB together with Acp using the selected 
dataset, while we correct the data for the detector- 
induced component using coefficients derived from a con- 
trol sample. 



IV. CORRECTION OF DETECTOR-RELATED 
ASYMMETRIES 



III. SIGNAL YIELD AND ASYMMETRY 
EXTRACTION 

For each mode the signal yield is extracted using a 
binned maximum likelihood (ML) fit to the distribution 
of the invariant mass m(K®h) for the selected D^s can- 
didates. The total probability density function (PDF) is 
the sum of signal and background components. The sig- 
nal PDF is modeled as a sum of two Gaussian functions 
for the D^s —> K^K^ 1 modes, and as a single Gaussian 

function for the Df — > K^tt^ mode. The background 
PDF is taken as the sum of two components: a distribu- 
tion describing the invariant mass of mis-reconstructed 
charm meson decays, and a combinatorial background 
modeling the mass distribution from other sources. For 
the D ± -> K°K^ (Df -> K ^) mode the charm back- 
ground is mainly from the tail of the invariant mass distri- 
bution for Df -> K°K± (D± -> K^) candidates. For 
the Df 

background originates mainly from D 
cays for which the ir^ 1 is misidentified as a K ± . As- 
signing the wrong mass to the pion shifts the recon- 
structed invariant mass, and the resulting distribution 
is a broad peak with mean value close to the Df mass. 
For each mode, the invariant mass distribution due to 
charm background is modeled using a histogram PDF ob- 
tained from a MC sample of simulated charm background 
decays. The combinatorial background is described by 
a first (second)-order polynomial for the Df — > 
mode (D ± ->■ and Df -> K^K* modes). The 

fits to the m(K° s h) distributions yield (159.4 ±0.8) x 10 3 
D ± -> K° S K^ decays, (288.2 ± 1.1) x 10 3 Df -> K^K^ 
decays, and (14.33 ± 0.31) x 10 3 Df -> decays. 
The data and the fit results are shown in Fig. [I] All of 
the PDF parameters are extracted from fits to the data. 

For each channel, we determine Acp by measuring the 
signal yield asymmetry A defined as: 



K[iK ± mode, the mis-reconstructed charm 

± - Kfr* de- 



.4 



N D+ -N D - 

(>) {•) 

N D+ +N D - 



(2) 



is the number of D^(D^) decays 

determined from the fit to the invariant mass distribu- 
tion. The asymmetry A contains two contributions in 
addition to Aqp, namely the forward-backward (FB) 



We use a data-driven method, described in detail in 
Ref. PP, to determine the charge asymmetry in track 
reconstruction as a function of the magnitude of the 
track momentum and its polar angle in the laboratory 
frame. The method exploits the fact that T(45) — > BB 
events provide a sample evenly populated with positive 
and negative tracks, free of any physics-induced asymme- 
tries. The off-resonance momentum distribution is sub- 
tracted from the on-resonance one, to remove any con- 
tribution from continuum, for which there is a FB asym- 
metry in the CM frame. This sample is used to com- 
pute the detector-related asymmetries in the reconstruc- 
tion of charged-particle tracks. Starting from a sample of 
50.6 fb" 1 of data collected at the T(4S) resonance and an 
off-resonance data sample of 44.8 fb" 1 , we obtain a large 
sample of charged-particle tracks, and apply the same 
charged pion or kaon track selection criteria used in the 
reconstruction of the D 



± 

(«) 



K^K* and Df -> K u s it 



modes. Then, after subtracting the off-resonance contri- 
bution from the on-resonance sample, we obtain a sample 
of more than 120 million pion candidates, and 40 million 
kaon candidates, originating from T(45) decays. We use 
the full off-resonance sample and an equivalent luminos- 
ity for the on-resonance sample, because, due to the sub- 
traction procedure, including additional data in the on- 
resonance sample does not improve the statistical error 
on the correction ratios mentioned below. These candi- 
dates are then used to compute the efficiency ratios for 
positive and negative pions and kaons. The ratio values 
and their statistical errors for pions and kaons are shown 
in Fig. [2] and Fig. [3J respectively. For the D^ s) -> K°K~ 

(D~ —> K®7r~) modes, the D7^ (D~) yields, in intervals 
of kaon (pion) momentum and cosine of its polar angle, 
cos 9, are weighted with the kaon (pion) efficiency ratios 
to correct for the detection efficiency differences between 
K + and K~ (n + and Momentum and cosine of 

its polar angle intervals are not uniform in order to have 
similar statistics, and therefore similar correction uncer- 
tainty, in each interval. Interval sizes vary from (0.05 
GeV/c, 0.06) to (4.4 GeV/c, 0.96), where the first num- 
ber is the momentum interval, and the second its cosine 
of polar angle interval. The largest correction is approx- 
imately 1% for pions and 2% for kaons. After correcting 
the data for the detector-induced component only A FB 
and Acp contribute to the measured asymmetry A. 
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FIG. 1: Invariant mass distribution for (a) D ± — > KsK ± , (b) Df — ¥ K^K^ 1 , and (c) Df — > K° s % ± candidates (points with 
error bars). The solid curve shows the result of the fit to the data. The dashed curve represents the sum of all background 
contributions, while the dotted curve indicates combinatorial background only. 
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FIG. 2: (top) The ratio between the detection efficiency for 
7r + and ir~ , and (bottom) the corresponding statistical errors. 
The values are computed using the numbers of tv + and n~ 
tracks in the selected control sample. 
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FIG. 3: (top) The ratio between the detection efficiency for 
K + and K~ , and (bottom) the corresponding statistical er- 
rors. The values are computed using the numbers of K + and 
K~ tracks in the selected control sample. 



V. EXTRACTION OF A C p AND A FB 

Neglecting higher-order terms that contain Acp and 
A-p-g, the resulting asymmetry can be expressed simply 
as the sum of the two. Given that A FB is an odd func- 
tion of cos 9* D , where 6* D is the polar angle of the Dj\ 
candidate momentum in the CM frame, Aqp and A FB 
can be written as a function of | cos6* D \ as follows: 

^(Icos^D^ ^ 1008 ^ 1 ^^ 1008 ^ (3) 
and 



where A(+|cos0£,|) (A(—\ cos 9* D \)) is the measured 
asymmetry for the candidates in a positive (neg- 
ative) cos 6^, interval. 

A simultaneous ML fit to the Dfs and D7-. invari- 

( s ) (s) 

ant mass distributions is carried out to extract the sig- 
nal yield asymmetry in each of ten equally spaced cos 9* D 
intervals, starting with interval 1 at [—1.0,-0.8]. The 
PDF model that describes the distribution in each sub- 
sample is the same as that used in the fit to the full 
sample, but the following parameters are allowed to float 
separately in each sub-sample (referred to as split pa- 
rameters): the yields for signal, charm background and 
combinatorial candidates; the asymmetries for signal and 
combinatorial candidates; the width, and the fraction of 
the Gaussian function with the larger contribution to the 
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signal PDF; and the first-order coefficient of the polyno- 
mial that models the combinatorial background. For the 
D ± — > K®K ± mode the yields for the charm background 
candidates in intervals 1, 2, and 3 were fixed to to ob- 
tain a fully convergent fit. Since interval 10 contains the 
smallest number of candidates, we use a single Gaussian 

function to model the signal PDF for the Df , — > K^K ± 

° ( s ) s 

modes. For the CP asymmetry of charm background 
candidates we use the same floating parameters as for 
the signal candidates, because the largest source of CP 
asymmetry for both samples is due to CPV in K° — K° 
mixing. For the Df — > Rgir mode, where the primary 
charm background channel, D ± — > R^tt^, has the same 
magnitude but opposite-sign asymmetry due to K°—K° 
mixing, we use a separate parameter for the asymmetry 
of the charm background candidates. To achieve a more 
stable fit, if the fit results for a split parameter are sta- 
tistically compatible between two or more sub-samples, 
the parameter is forced to have the same floating value 
among those sub-samples only. For the Df —> R^Tt 
mode the width of the first Gaussian function for the sig- 
nal PDF is set to the same floating value in intervals 1,2, 
3, and 4. The first-order coefficient of the polynomial de- 
scribing the combinatorial background is set to the same 
floating value in intervals 4 to 8 (D ± —> K I ^K ± ), in in- 
tervals 4 to 8 (Df K^K^), and in intervals 2 to 7 
(Df -> K^). The final fit contains 70, 80, and 64 free 
parameters for the D ± -> K° S K ± , Df -> R Q S R ± , and 
Df — > K®ir modes, respectively. 



The Acp and A FB values for the five | cos 6* D | bins are 
shown in Fig. [4] for the three decay modes. The weighted 
average of the five Aqp values is (0.16 ± 0.36)% for the 
D ± -> R°R ± mode, (0.00 ± 0.23)% for Df -+ R°R ± , 
and (0.6 ± 2.0)% for Df — > R^tt^, where the errors are 
statistical only. 



We perform two tests to validate the analysis proce- 
dure for each channel. The first involves generating 5000 
toy MC experiments with a statistics equal to data us- 
ing the PDF and the parameters obtained from the fit 
to data. After extracting Aqp from each experiment, for 
the D ± -¥ K^K ± and Df -> K°K ± modes, we deduce 
from the mean of the Aqp pull distributions the presence 
of a small bias in the fitted value of each fit parameter 
(the means are -0.036±0.014 and +0.041±0.014, respec- 
tively) . To account for this effect we apply a correction to 
the final values equal to +0.013% for the D ± -> K^K* 1 
mode, and -0.01% for the Df -> K°K ± mode. The A C p 
pull distributions show that the fit provides an accurate 
estimate of the statistical error for all the modes. The 
second test involves fitting a large number of MC events 
from the full BABAR detector simulation. We measure 
Aqp from this MC sample to be consistent with the gen- 
erated value of zero. 



VI. SYSTEM ATICS 

The main sources of systematic uncertainty are listed 
in Table [TT] for each decay mode, together with the overall 
uncertainties. The primary sources of systematic uncer- 
tainty are the detection efficiency ratios used to weight 
the .D7j yields, and the contributions from mis-identified 
particles in the data control sample used to determine the 
charge asymmetry in track reconstruction efficiency. 

The technique used to remove the charge asymme- 
try due to detector-induced effects produces a small sys- 
tematic uncertainty in the measurement of Aqp due to 
the statistical error in the relative efficiency estimation. 
This systematic uncertainty depends only on the type of 
charged particle (pion or kaon) in the final state, and 
not on the initial state. To estimate the systematic un- 
certainty on Acp resulting from this source, the relative 
charged-particle efficiency in each interval of momentum 
and cos 6 is randomly drawn from a Gaussian distribu- 
tion whose mean is the nominal relative efficiency in that 
interval, and where the root-mean-squared (r.m.s.) de- 
viation is the corresponding statistical error. For each 
mode, we generate 500 such charged-particle relative- 
efficiency distributions, and use them to obtain 500 Aqp 
values, following the procedure described earlier to deter- 
mine the nominal value of Aqp- The r.m.s. deviation of 
these 500 values from the nominal Acp is taken to be the 
systematic uncertainty. For the D^ — > K®K ± modes, 

the estimated systematic uncertainty is 0.23%. For the 
Df -> K^ir* mode, we assign the same systematic un- 
certainty, 0.06%, as that estimated for the D ± — > R^tt 
mode in Ref. [1]. 

The small fraction of misidcntificd particles in the 
generic track sample can introduce small biases in the es- 
timation of the efficiencies, and subsequently in the Aqp 
measurements. Because of the good agreement between 
data and MC samples, we can use the simulated MC can- 
didates to measure the shift in the Aqp value from the 
fit when the corrections are applied, and when they are 
not. Again, this contribution depends only on the type of 
the charged-particle track. Hence, for the Df — > 
mode, we assume the same shift obtained in Ref. [T], 
namely +0.05%. By fitting the Df -> K° S R ± MC sam- 
ple when the corrections are applied, and again when not, 
we obtain a shift of +0.05% and we assume this for both 
the Df -> K°K± and D ± -> K°K ± modes. For all 
the modes, we shift the measured Aqp by this correction 
value and then, conservatively, include the magnitude of 
this shift as a contribution to the systematic uncertainty. 

Using MC simulation, we evaluate an additional sys- 
tematic uncertainty of ±0.01% due to a possible charge 
asymmetry present in the control sample before applying 
the selection criteria. Another source of systematic un- 
certainty is due to the choice of the cos 9* D interval-size 
in the simultaneous ML fit. The systematic uncertainty 
is taken to be the largest absolute difference between the 
nominal Aqp extracted using ten cos 9* D intervals and 
that obtained when the fit is performed using either 8 
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FIG. 4: CP asymmetry, A C p, for (a) D i 



K^K ± , (b) Df IsTgiir*, and (c) £>± Kgvr* as a function of | cos6lo| in the 
data sample. The solid line represents the central value of Acp and the gray band is the ±1 a interval, both obtained from a 
X 2 -minimization assuming no dependence on | cos#J>|. The corresponding forward-backward asymmetries, A FB , are shown in 
(d), (e), and (f). 

TABLE II: Summary of the systematic uncertainty contributions for the Acp measurement in each mode. The values are 
absolute uncertainties, even though given as percentages. The total value corresponds to the sum in quadrature of the individual 
contributions. 



Systematic uncertainty 


D ± -> K° S K ± 


Df -> Ai^± 


Df -> AItt* 


Efficiency of PID selectors 


0.05% 


0.05% 


0.05% 


Statistics of the control sample 


0.23% 


0.23% 


0.06% 


Misidentified tracks in the control sample 


0.01% 


0.01% 


0.01% 


cos 9* D interval size 


0.04% 


0.02% 


0.27% 


K°—K° regeneration 


0.05% 


0.05% 


0.06% 


Kg — interference 


0.015% 


0.014% 


0.008% 


Total 


0.25% 


0.24% 


0.29% 



or 12 intervals in cos9* D . This is the dominant source 
of systematic uncertainty for the Df — > AT^Tr* mode, as 
shown in Table HT1 

We also consider a possible systematic uncertainty due 
to the regeneration of neutral kaons in the material of the 
detector. The K° and K° mesons produced in the decay 
processes can interact with the material in the tracking 
volume before they decay. Following a method similar to 
that described in Ref. [IB], we compute the probability 



for a K° or a K° meson to interact inside the BABAR 
tracking system, and estimate systematic uncertainties 
of 0.05% (Df s) -> A^A:*) and 0.06% (Df -> A^vr*). 

Although the intermediate state is labelled as a AT°, 
we apply a correction term to the measured Acp to in- 
clude the effect of Kg— K® interference in the intermediate 
state [19] . This correction term depends on the proper 
time range over which decay distributions are integrated, 
and on the efficiency of the reconstruction of the 7r + 7r~ 
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TABLE III: Summary of the Acp measurements. Where reported, the first uncertainty is statistical, and the second is 
systematic. 





d ± -+ kIk* 






Acp value from the fit 


(+0.155 ±0.360)% 


(0.00 + 0.23)% 


(+0.6 + 2.0)% 


Correction for the bias from toy MC experiments 
Correction for the bias in the PID selectors 
Correction for the Kg — K\ interference (AAcp) 


+0.013% 
-0.05% 
+0.015% 


-0.01% 
-0.05% 
+0.014% 


-0.05% 
-0.008% 


Acp final value 


(+0.13 + 0.36 + 0.25)% 


(-0.05 + 0.23 + 0.24)% 


(+0.6 + 2.0 + 0.3)% 


Acp contribution from K° — K° mixing 


(-0.332 + 0.006)% 


(-0.332 + 0.006)% 


(+0.332 ± 0.006)% 


Acp final value (charm only) 


(+0.46 + 0.36 + 0.25)% 


(+0.28 + 0.23 + 0.24)% 


(+0.3 + 2.0 + 0.3)% 



final state as a function of proper time. We compute 
the reconstruction efficiency distribution as a function of 
proper time using MC truth-matched K® decays after 
the full selection. Following the method in Ref. [19] we 
estimate the asymmetry-correction term AAqp defined 
as: 



atic. In the last row of the table, we also report the Aqp 
values after subtracting the expected Acp contribution 
for each mode due to K°~K° mixing. The results are 
consistent with zero, and with the SM prediction, within 
one standard deviation. 



AA 



CP 



a corr 



A 



fit 
CP' 



(5) 



where A^p is the value obtained from the fit and A^p is 
the corrected value. The correction terms are reported in 
Table [Il*l*| and, to be conservative, we include their abso- 
lute values as contributions to the systematic uncertainty 
estimates. We also estimate the correction factor for the 
D ± -> K° 7 r ± mode using the K® reconstruction effi- 
ciency distribution after the selection detailed in Ref. pQ , 
and obtain the value +0.002%. All these corrections are 
rather small, even compared to those estimated in a sim- 
ilar analysis [2D]. The smaller values of the corrections 
in the present analysis are due to the improved efficiency 
for Kg mesons with short decay times that we obtain by 
applying the requirement on the decay length divided by 
its uncertainty, rather than on the decay length alone. 



VII. CONCLUSION 

In conclusion, we measure the direct CP asymmetry, 
A CP , in the D ± -> K°K±, Df -> K^K*, and Df -> 
K^ir^ modes using approximately 159 000, 288 000, and 
14 000 signal candidates, respectively. The measured 



Aqp value for each mode is reported in Table III where 



the first errors are statistical and the second are system- 
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